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Abstract. We propose a measurement-based method to produce a maximally- 
entangled state from a partially-entangled pure state. Our goal can be thought of as 
entanglement distillation from a single copy of a partially-entangled state. The present 
approach involves local two-outcome weak measurements. We show that application of 
these local weak measurements leads to a probabilistic amplification of entanglement. 
In addition, we examine how the probability to find the maximally-entangled state 
is related to the entanglement of the input state. We also study the application of 
our method to a mixed initial state. We show that the protocol is successful if the 
separable part of the mixed initial state fulfils certain conditions. 
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1. Introduction 

Measurements on quantum systems as well as the coherent manipulation via unitary 
operations are key ingredients to the implementation of interesting quantum state- 
engineering protocols. An essential feature of quantum measurements jTJ [2j [3], |4] is that 
these significantly affect the static and dynamic properties of the system. This allows 
phenomena not predicted by classical mechanics, such as the quantum Zeno effect [51 [6] . 

Many measurement-based schemes have been studied including: one-way quantum 
computation [7], protection or recovery of a quantum state in a noisy channel [H 
El Uni EE], preparation of an entangled state [12], and measurement-based quantum 
control [T3im] . Among these approaches, the use of weak (or general) measurements has 
recently drawn considerable attention. General measurements are the generalization of 
von Neumann measurements and are associated with a positive-operator valued measure 
(POVM) [U [2j [3]. Their influence on a quantum system is, in general, more moderate 
than von Neumann measurements. Thus, in some cases, the amount of information 
extracted using weak measurements can be tunable. This feature is useful to implement 
quite interesting quantum protocols, such as reversing measurements [15j [16] and weak- 
value measurements [HI HE] . 

In this paper, we study a method to prepare a maximally-entangled state from 
a pure initial state using weak measurements. Our goal can be considered to be 
entanglement distillation with a single copy of a partially-entangled state. The present 
approach involves local two-outcome weak measurements. We show that the application 
of local weak measurements leads to a probabilistic amplification of the entanglement. 
Furthermore, we examine the relationship of our scheme with an entanglement- 
concentration protocol [HI ED]- The present study shows several interesting roles of 
a weak measurement. 

The paper is organized as follows. First, we explain our setting in section El Section 
[3] is the main part of this paper. We propose a method to produce a maximally- 
entangled state from a partially-entangled pure state with the application of local weak 
measurements. The key idea is to tune the measurement strength at each measurement 
step. Our approach developed in section [3] assumes that the input state is a pure state. 
In order to consider more realistic situations, we examine the case when the input state 
is a mixed state in section HI Then, we find that the method for a pure state still leads to 
entanglement amplification if the separable part of the initial mixed state fulfills certain 
conditions. Furthermore, we discuss the extendability of this approach in section [5j 
Section [6] is devoted to a summary of the results. 

2. Setting 

2.1. Linear entropy 

Let us consider a bipartite system composed of a two- and a cHevel system (d > 2). 
The two-level system is called system A, while the ci-level one is called system B. We 
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will apply some operations only to system A and never touch system B. Since system 
A is a single-qubit, the treatment of control processes is rather simple. We have a pure 
state in this composite system, 



with a 2 + f3 2 = 1, a > 0, (3 > 0, = 5 ij} and (4>i\(f)j} = d~ij (i,j = 0,1). We also 
note that a and (3 are real parameters. Mathematically, we can always find this form 
of a vector with the Schmidt decomposition |3J. The quantum correlation of is 
characterized by the purity of the reduced density matrix pa = Ttb(|'0)(V , I)) as seen i n > 
e.g., [4J. We consider the linear entropy 



When S(ip) = 0, we find that either a = or (3 = 0. Hence, there is no correlation 
between the two subsystems. In contrast, when S(ip) > 0, has a quantum correlation, 
i.e., entanglement. In particular, S(ip) takes its maximum value 1 when pa is a 
maximally- mixed state (i.e., when a = (3). This indicates that we have a maximally- 
entangled state. Thus, the linear entropy S(ip) can be regarded as an entanglement 
measure. 

2.2. Motivation 

Entanglement is one of the key resources of quantum information processing such as 
quantum teleportation, quantum dense coding, quantum key distribution, and so on |21j. 
To develop a way to create and protect maximally-entangled states is quite important 
because most quantum protocols rely on their existence. Thus, an interesting physical 
issue is that ([T]) is a partially-entangled state (i.e., a ^ (3 and a, (3 > 0) and one 
might desire to amplify its entanglement. Typically, one can encounter such an issue 
when one tries to create a maximally-entangled state in the presence of experimental 
imperfections. Our aim is to find an efficient and practical way to transform \ip) into 
a maximally-entangled state. The most straightforward way might be to use controlled 
unitary operators between the two subsystems. Entanglement purification with many 
copies of \if)) is also a candidate to accomplish this goal, as seen, e.g., in [21]. However, 
these methods can be costly, e.g., either requiring properly controlled gates, or preparing 
a sufficient number of copies, or two-body measurements. Here we wish to propose a 
simpler method to achieve the same goal. Namely, we want to create a maximally- 
entangled state from a partially-entangled state using a simpler experimental setting. 
Thus, we consider the case when one can use only local operations and a single copy of 



We now illustrate our task more explicitly. Without loss of generality, we assume 
that (3 > a in ([1]). Let us consider a completely positive and trace-preserving (CPTP) 
map [2 J 




(1) 



^) = 2[l-Tr A (p A )]=4a 2 /3 2 . 



(2) 




(3) 
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with a set of linear operator on C 2 , Ke = <8> Id and "^2 1 K\Ki — I 2 <S> Id- The 
identity operator on C 2 (C d ) is denoted by I 2 (Id)- The entanglement of p out must 
be smaller than the entanglement of p m because this map is composed of only local 
operators. However, a portion of p out , e.g., Ki\ip)(ifj\Kj /Tr (Ki\ip) (ip\Kj) can be a more 
strongly entangled state than p m [21] . The local map ([3]) allows entanglement distillation 
with probability Tr (Ki\tfj)(ip\Kj). Thus, the problem reduces to seeking a set of K^s 
and how to apply them. 

Our goal can be thought of as entanglement distillation from a single copy 
of a partially-entangled state. Bennett et al. [19] suggested a scheme to purify a 
maximally-entangled state from a partially-entangled pure state using local partial- 
collapse measurements [15] , i.e., = v /p|l)(l| and K 2 — |0)(0| + y/1 — p\l) (1|, with 
< p < 1. Kwiat et a/. [20] reported the experimental realization of this scheme. 
Although one has to discard the quantum state when the outcome is K^K±, one 
obtains a maximally-entangled state when the outcome is K^K^- Several other 
entanglement distillation procedures have been proposed (e.g., [22]) but using many 
copies of a partially-entangled state. The demonstration of the methods with two 
copies of partially-entangled states was reported in, e.g., [231 121] • The basic setting 
in this paper is the same as in JT9J [20] (i.e., the use of a single copy of a pure state and 
local measurements). Thus, hereafter, we will use the term "distillation" for describing 
our protocol although we do not treat a large ensemble of mixed states but just a single 
copy of a partially-entangled pure state. 

We propose a method for entanglement amplification with local weak measurements 
and a single copy of a partially-entangled pure state. If the measurement is successful, 
we obtain a maximally-entangled state. We stress that our approach requires a small 
cost because a lot of copies and two-body measurements are not needed. The type of 
the weak measurements used in our protocol is different from the one in the proposal 
by Bennett et al. (so-called Procrutean method) [19] . and was demonstrated in, e.g, 
[17] . Although our approach is similar to the Procrutean method, the present method 
does not completely destroy entanglement even if the measurement fails. Thus, we can 
repeatedly apply weak measurements until we obtain a maximally-entangled state. We 
also remark that when one has a physical system in which the implementation of the 
weak measurements in [19] is difficult, one can use the present approach as an alternative 
method for entanglement distillation of a single copy of a pure state. 

3. Entanglement amplification with local weak measurements 

Now we show a method to distill entanglement from a partially-entangled pure initial 
state. The basic part of the scheme is the weak measurements on system A given by 

M ± (6) = y^|0)(0| + y^p|l>(l|, (4) 

m|m+ + \r m = j 2 , (5) 
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Figure 1. Schematic diagram of our approach to distill a maximally-entangled state 
from a partially-entangled initial state = a\0)\(f>o) + The measurements 

are performed only on system A. Each box indicates the measurement apparatus with 
the measurement strength e„ (n = 1, 2, . . . , N), see (U]). When we obtain the outcome 
corresponding to M±M+, we find a maximally-entangled state. 




Figure 2. Application of our approach to a|O)|0o) + with a 2 = 0.4 and 

P 2 = 0.6. The horizontal axis is the number of measurements performed. The 
probability to obtain a maximally-entangled state is shown in (a). The measurement 
strength is shown in (b). The sum of the net probabilities, which is defined as (fl"2"]) . is 
shown in (c). The dotted line indicates the total success probability, 2a 2 = 0.8. See 
(Tl3l) and (|T7| . A higher total success probability could be easily obtained by choosing 
2a 2 close to one. 
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with measurement strength e (0 < e < 1). The set {M^M a } a=+ _ is a POVM [IJ [2J E] 
on C 2 . These measurements with tunable e can be implemented in various physical 
systems, as seen, e.g., in [I7J [25j EH] • We substitute M+(e) and M_(e), respectively, into 
K± and K% in The resultant state is given by p out = p + \ip + ) \ +p-\ip-)(il ! - \ with 
|-0±) = {M± <g) Id)\i>) / y/P± and p± = (^|(Mj_M± ® The important point here 

is that 1-0+) is a maximally-entangled state when a 2 (l + e)/2 = /3 2 (1 — e)/2. Assuming 
that (3 > a, we find that the parameter e for obtaining a maximally-entangled state 
becomes 

e = p 2 -a 2 = y/l-S(i>). (6) 

In contrast to a partial-collapse measurement [15], the present measurement operators 
do not completely destroy the entanglement of the input state even if the measurement 
fail. We also note that both the present approach and the Procrustean method [19] 
require a priori information of the initial state (i.e., a and f3). 

On the basis of the above arguments, we propose a probabilistic method to make 
a maximally-entangled state from with the successive application of local weak 
measurements, as shown in figure [TJ We prepare the N weak measurement apparatus 
described by (jlj) and fl5]). We stop the measurement process once we obtain the outcome 
M|M + . Otherwise, we move to the subsequent measurement. The input state in the 
nth weak measurement apparatus is 

\^n) = r J M-(e w _i) ® IdMn-l), (7) 

P-(en-i) = (V'„-i|M_(e n _ 1 ) t M_(e n _ 1 ) <g> I d \^i), (8) 

if n > 2 and |?/>i) = with > a. The expression of the weak measurements (jlj) 
allows the form invariance of \ip n ), i.e., \ip n ) = a n |O)|0 o ) + /3 n ,|l)|0x) where a 2 + /3 2 = 1 
and a n , f3 n > 0. In addition, we find that /3 n > a n if /3 n _i > and e n > 0. Then, 
when the measurement strength obeys the recurrence relation 

en = , *T X , 2 (n>2), c^V^W, (9) 

the state corresponding to the outcome M^(e n )M + (e n ) ® J d is the maximally-entangled 
state (|O)|0o) + |l)|0i))/v2- Indeed, the solution of (J9j) satisfies jB}, i.e., e n = 
y/l — S(ip n ). The probability to obtain the outcome M\_M + at the nth measurement 
apparatus is 

Pn = III- (e n ) 2 ]. (10) 

When /3 < a, we can obtain the same result replacing the role of M + with M_. Since 
e n obeys the recurrence relation (Q, the information required is only S(ijj). 

Figure |2] shows the probability P n and the measurement strength e n for a specific 
input state. After a few steps, the probability decreases and approaches 0, and the 
measurement strength increases and approaches 1. Actually, the recurrence relation 
(EJ) indicates that e n > e n _i, and linin^ooen = 1. Therefore, the weak measurements 
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asymptotically approach the von Neumann measurements |0)(0| and |1)(1| for large n. 
Then, the state corresponding to the outcome M [ _M_ ® becomes 



Figure 3. 

represents 
defined as 
He). 



Total success probability versus initial entanglement. The horizontal axis 
the linear entropy S(\tp)) of the state The total success probability is 
the asymptotic value of the sum of the net probability, as shown in figure 




Initial entanglement 



Let us now characterize better the present approach. The nth measurement is used 
when the previous weak measurement apparatus shows the outcome M [ _M_. Thus, we 
may define the net probability to find a maximally entangled state at the nth step as 

n-l 

PT = Pn\[{l-Pk), (11) 

fc=l 

if n > 2 and Pi ct = P x . We find that P£-i > P™ 1 - The sum of the net probabilities is 
a characteristic quantity of this scheme. We define P^ as 



Epnet 



(12) 



m=l 



n-l 



+ 



;i - P n )Pn-i if n > 2 and 



This quantity obeys the recurrence relation P^ 
PI = P\. Figure Ws) shows that P^ approaches a specific value for large n. We confirm 
that this asymptotic value is 2a 2 by numerically calculating P^ for the various initial 
states. We call this asymptotic value the total success probability of the approach. 
Figure [3] shows that the total success probability is large when a is close of \ [or S^VO is 
close to its maximum value]. Here, we show an alternative evaluation of the total success 
probability. We consider an asymptotic map corresponding to the entire protocol. First, 
we focus on the fact that the entire protocol in the limit n — > oo can be regarded as a 
CPTP map 

|V>)<VI ^ P* = WmesX^mesI + (1 - P*)\X){X\, (13) 

where |^mes) = (|O)|0o) + |l)|0i))/\/2, and |x) = |l)|0i). The positive coefficient 
corresponds to the total success probability. This map must be a composite map of the 
weak measurements (j3]). Therefore, we may write p* as 

p* = (M, t0 ® I d )|V>)(^|(M* )CT ® h)\ (14) 



M, 



*.± 



l±6c. 



|0)(0| + 



1±5 1 



1)(1| 



(15) 
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with \5q\, I 5x| < 1. The measurement operators M*,± correspond to generalized weak 
(partial) measurements [26 [ 127] . We find that (fT4"j) represents a convex combination 
between |^mes)(^mes| and \x)(x\ if 

(M*,+ ® I d )\ip) cx I^mes), (M» _ ® J d )|-0> oc |x). (16) 
Let us put 5 = 1 and fly/ (1 + Si)/2 = a. Then, we find that (M*,+ <8> IcdO*) — 
V^al^MEs)- Thus, we show that 

P* = 2a 2 . (17) 

In other words, we decomposed the map (TT51) with PJ, as pointed out in [30] • The 
measurement operators M*,+ and M*,_ for 5q = 1 was experimentally demonstrated as 
a partial-collapse measurement in [TBI 1201 EH]- We remark that the asymptotic map 
can be regarded as a single-shot partial measurement and is equal to the Procrustean 
method in [19J. Actually, the total success probability P* is the same as the yield of 
maximally-entangled states in [19]. 

Let us compare the present entanglement amplification protocol to the Procrustean 
method in [19] from the viewpoint of physical implementation. The present approach 
requires the implementation of weak measurements that are symmetric or balanced with 
respect to the measurement strength, as seen in fll]). Such symmetric weak measurements 
can be implemented in various physical systems including linear optical and solid- 
state qubits, as shown, e.g., in [29]. In contrast to the symmetric measurements, the 
Procrustean method uses non-symmetric ones. Namely, one uses ( IT5l) with 5q = 1 and 
Si 7^ —1. Although the non-symmetric partial measurements are naturally realized in, 
e.g., Josephson phase qubits [28] . they are not always available in a general physical 
system. Thus, when one has a physical system in which the implementation of the type 
of non-symmetric partial measurements that are used in [19] is difficult, one can use the 
present approach as an alternative method for entanglement distillation of a single copy 
of a pure state. 

4. Mixed states 

We now extend our approach to a mixed initial state. Hereafter, we consider a two-qubit 
system (i.e., d = 2). The study of a mixed initial state is motivated by experimental 
considerations, as well as theoretical interest. One may desire to amplify entanglement 
of a quantum state in the presence of decoherence, for example. We remark that one 
never distills a maximally-entangled state from a single copy of a mixed state with 
local purification protocols, as shown by Linden et al. [31] and Kent [32] . Hence, our 
goal in this section is to seek a protocol of increasing the entanglement of an output 
state compared to a mixed initial state, not purifying a maximally-entangled state. This 
section mainly focuses on a single-shot measurement protocol, not a repeated application 
of weak measurements. 

The initial pure state ([1]) is replaced with the mixed state 

p = \p B + (l-\)\iP)(iP\, (18) 
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where = a|00) + and < A < 1. This is a convex combination between a 

separable state p s and a partially-entangled pure state \ip). The coefficients a and 
are real and positive numbers in the same way as (CD). The weight of the separable 
part is quantified by A. We note that an arbitrary density matrix on C 2 <E> C 2 can 
be written by the form ( fl8l) . and this decomposition is uniquely determined [35] . We 
use the concurrence [331 El] to evaluate the entanglement of a mixed state. Hereafter, 
the concurrence of p is written as C(p). We also use another quantity to estimate 
entanglement in addition to the concurrence. When a mixed state is given as the form 
(|T8|) . its entanglement can be characterized by its pure-state part (1 — X)\xjj){^\ [35J. 
Thus, an intuitive "measure" of entanglement is defined as E(p) = (1 — \)S(ip), with 
the linear entropy ([2]). Using both C(p) and E(p), we examine the entanglement 
amplification protocol for (Tl8l) . Hereafter, we assume that A 7^ and (3 > a. 

We use the same local two-outcome weak measurement in the case of a pure initial 
state. As seen in section [3j performing a single measurement whose measurement 
operator is given by (j4j), we obtain p out = + R-p~, with R± = Xw St ± + (1 — A)p±, 

p± = (M ± <g) J 2 )p(M| (g) I 2 )/R±, and w Sj ± = tr [(M£M ± (g) i^As]- O ur a i m * s ^° increase 
the entanglement of p + , compared with p. This density matrix has the same form as 



P+ 



^°i>)»M°h) + {l _ K)m)l M, (19) 



with 



A + = ^±A. (20) 



We examine E(p) and E(p + ) to obtain a simple criterion for the entanglement 
amplification. Since p s is a separable density matrix and the weak measurement is local, 
the first term in ( fl9i) is the separable part of p + . Therefore, E(p + ) = (1 — \ + )S(i[) + ). As 
shown in section [3], 1^+) is a maximally-entangled state when the measurement strength 
satisfies (EI). It means that S(ip + ) > S(ip). To ensure E(ip + ) > E(ip), we require A + < A. 
This relation leads to 



A, z < -y/l-S(4>) = a 2 - /3 2 , (21) 

where we have used ([6]). We have defined A S)Z as A s<z = tr [{a z <8> l2)p s ], where 
o~ z = |0)(0| — Therefore, we expect that the method for a pure initial state 

works even for a mixed initial state if this inequality is satisfied. We will check this 
simple criterion for specific examples using calculations of the concurrence. Then, our 
numerical calculations of the concurrence will show that ( 12T1) is the sufficient condition 
for C(p+) > C(p). 

4-1. Example 1: Pure dephasing 

Let us consider the case when one of the subsystems (e.g., system A) travels via a noisy 
channel after the preparation of \ip), for example. Assuming that the system undergoes 
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pure dephasing, the corresponding CPTP map is written as p h-> (1 — u)p + u[a z ® 
I%)p(pz <8> I2), with < u < 1/2. Thus, the initial mixed state becomes 

p( pd ) = « 2 |00)(00| + u|00)(ll| +/i|ll)(00| + /3 2 |11)(11|, (22) 

with // = a/3(l - 2m). We can rewrite (J225 as p( pd > = [1 - (u/a/3)](a 2 |00)(00| + 
/9 2 |11)(11|) + (p/ce/3)\il))(if}\. Non-classical features of this sort of density matrices 
are characterized well by the off-diagonal element /^[36]. In fact, the concurrence 
is C(p( pd )) = 2/i. Let us check the inequality ( 12T1) . Since A BiZ = a 2 — (3 2 and 
y/l - S(ip) = f3 2 -a 2 , the left hand side of is zero and the inequality is satisfied. 
Thus, we expect that the approach for a pure initial state leads to an amplification of 
the entanglement of p( pd \ We calculate the concurrence to confirm that this prediction 
is correct. The concurrence of p+ d ^ oc (M+ <g> J 2 )p^ pd ^(M^ Cg> J 2 ) is 

C(p pd , + ) = -^->2p = C(p^). (23) 

Therefore, entanglement amplification is achieved. Furthermore, we can find that the 
repeated application still works. The input density matrix for the nth measurement 
apparatus is the state corresponding to the outcome Ml(e n _i)M_(e n _i). We denote 
this state as p n . We note that pi = p^ pd \ When we obtain the outcome M|(e n )M + (e n ), 
the measurement is considered to be successful. We write the resultant state as pi . 
The measurement strength obeys the recurrence relation with e\ = v/l — S(ip). The 
application to p is straightforward because p n and pn have the same form as p (i.e., 
only the coefficients change). We can show that C{p < t > ) > C(p n ) and C(p n ) < C(p„_i). 
Eventually, we find that C(pi^) = C(p t ^_ 1 ). Therefore, at every value for the outcome 
we obtain a mixed state with the concurrence C(p^ ) = p/a(3. Remarkably, we 
find that u/a/3 > C(p( pd ^) because 2a(3 < 1. The resultant entanglement is bounded by 
p/a(3, which is regarded as a measure of the purity of p( pd \ 



4-2. Example 2: Amplitude damping 

The second example is associated with amplitude-damping. If this effect occurs on 
system A, the corresponding CPTP map is written as p ^ J2e=i 2(^1 <S> I 2 )p{E\ <S> I2), 
with E x = |0)(0| + Vu~\l)(l\, E 2 = y/T=u\0)(l\, and < u < 1. Thus, the mixed state 
becomes 

pM = (1 _ M )/3 2 |01)(01| + [1 - (1 - u)p 2 }\ij)(ij\, (24) 

with = (a|00) + ^/uP\U))/y/l - (l-u)/3 2 . Assuming that 01/3 > a, we check the 

inequality flUD- Since yjl- S(ip) = (u(3 2 - a 2 )/(u(3 2 + a 2 ) > and A^ z = 1, fl2D is 
not satisfied. Thus, we expect that the method for a pure initial state does not work. 
A numerical calculation of the concurrence gives more quantitative characterization of 
our protocol, as seen in figure HI Varying the weight of the separable part (1 — u)(3 2 
and the linear entropy of the reduced density matrix of the pure-state part S(ip), the 
concurrences of p^ [left panel of figure BJa)] and p+ d ^ oc (M + ® J 2 )p ( - ad ' , ( j ^+ ® h) 
[right panel of figure IHa)] are calculated. Furthermore, the sign of their difference ( 
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Figure 4. (Color online) Characteristics of the entanglement amplification method 
for the mixed state associated with amplitude damping, (|24|) . The horizontal axis is 
the linear entropy of the reduced density matrix of the pure-state part, S(ip)- The 
vertical axis is the weight of the separable part, (1 — u)/3 2 . (a) Concurrence of the 
input state p( ad ) (left panel) and the output state pS* 1 (right panel). Using (j6]), the 

measurement strength is given by yJl-S$). (b) Profile of sgn[C(p ( + d) ) - C(p^)]. 
When this quantity is positive, the entanglement amplification is achieved (red-colored 
region). Otherwise, the protocol fails (blue-colored region). 




0.3 0.6 0.9 0.3 0.6 0.9 0.3 0.6 0.9 

Linear entropy S(ip) 



sgn[C(p+ ) — C(p( ad ))] ) is evaluated, as seen in figureH^b). Figure@Jb) shows that the 
method for a pure state leads to the reduction of the entanglement in wide parameter 
region. Nevertheless, when the purity of p^ is large, the entanglement is amplified. 

4-3. Example 3: Maximally-mixed state 

The parameter region in which the entanglement amplification for mixed states works 
depends on the value of A S)Z . Let us examine the case when the separable state in ffT8]) 
is the maximally-mixed state, for example. The initial mixed state becomes 

p(- d ) = A~I 2 ®I 2 + (l-A)|V>)<# (25) 

We find that A S)Z = 0. Therefore, the inequality ( J2TT) is not satisfied except for 
S(if)) = 1. In order to evaluate quantitative behaviors of the protocol, let us examine 
the concurrence. Figure G2a) shows the concurrences of p( rnd ) (left panel) and p^ nd ' ) oc 
(M + ®I 2 )p rd (Ml®I 2 ) (right panel). The sign of their difference (sgn[C(p5: nd) )-C(P (md) )] 
) is shown in figure \Mjo)- Again, we find that the method for a pure state fails in wide 
parameter region. However, compared with figure IHb), the entanglement amplification 
is achieved in much wider region. Thus, smaller values of A SjZ lead to higher success of 
the protocol for various parameters. We also find that from ( |2TI) negative values of A 8jZ 
are preferable for the entanglement amplification. 
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Figure 5. (Color online) Characteristics of the entanglement amplification method for 
the initial state (|25l) . The horizontal axis is the linear entropy of the reduced density 
matrix of a pure-state part, S(ip). The vertical axis is the weight of the separable 
part, A. (a) Concurrence of the input state p( rnd ) (left panel) and the output state 
^(r nd ) (j-jght panel). Using ((6]), the measurement strength is given by \/l — S(ip). (b) 
Profile of sgn[C(p^ nd ' ) ) — C( / o( md - ) )]. When this quantity is positive, the entanglement 
amplification is achieved (red-colored region). Otherwise, the protocol fails (blue- 
colored region). In the white-colored region, the concurrence does not change because 
the initial mixed state (|25|) is separable. 
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4-4- Systematic examinations 

For a more systematic analysis for mixed initial states, we randomly generate separable 
density matrices. Thus, we examine entanglement amplification for various density 
matrices given by (|T8l) . with fixed A separable state on C 2 ® C 2 is parametrized by 
15 real parameters. Using the Pauli matrices, we can express it as 

Ps = 4 ( h ® h + ^2 As ' i(Ti + Bs >^ 2 ® a i + C ^i a i ® a J ) ' ( 26 ) 

V i j i,j / 

with i,j = x,y,z and A Sti , B s j, C S)i j G M. The positivity and the separability of (|26|) 
give constraints with respect to the real parameters A Sji , B s j, and C Sji j G R The 
former is straightforwardly checked using Newton's formula (see, e.g., [37]). The latter 
is checked by the positive partial transpose criterion |38[ [39] . We calculate the mean 
values of the concurrences of p and p + with respect to randomly-generated separable 
states. Then, we evaluate sgn[C(p + ) — C(p)], where C represents the mean value 
of the concurrence. To examine the A s z dependence of the protocol, we randomly 
generate p s fixing A s as a specific value. We generate 10000 samples for each value 
of A S)Z using the Mersenne Twister method [40J. Figure [6] shows sgn[C(p + ) — C(p)], 
varying A SjZ . When this difference is positive (red-colored regions), the entanglement 
amplification is achieved. We note that in the white-colored region the entanglement 
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Figure 6. (Color online) Profiles of concurrence differences between the output and 
the initial states, varying A SyZ . Each point is a mean value over 10000 separable states. 
The horizontal axis is the linear entropy of the reduced density matrix of the pure- 
state part. The vertical axis is the weight of the separable part. See (fl8|) . When 
the differences are positive (red-colored regions), the entanglement amplification is 
achieved. Otherwise, the protocol fails (blue-colored regions). In the white-colored 
regions, the concurrence does not change because the mixed initial state is separable. 
In the right regions of the the vertical (green) lines, the inequality (|2Tj) is satisfied. 
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is not changed because the mixed initial states are separable. The vertical (green) 
lines (for A s>z = —0.95, —0.75, —0.55, —0.35) corresponds to the equality of f l2Tj) . 
S(i/j) = l-(A SyZ f. In the right regions of the lines, the inequality ( |2T|) is satisfied. Thus, 
figure [6] shows that ff2T]) is the sufficient condition for C(p + ) > C(p). Furthermore, when 
A SjZ is a large negative value, the entanglement is enhanced in wide parameter region. 

Summarizing the above arguments, the inequality fT2~Tj) is the sufficient condition for 
the entanglement amplification. In addition, the method for a pure intial state can still 
lead to entanglement amplification even for a mixed initial state if the separable part of 
the mixed state fulfills A SiZ < 0. The reason why negative values of A BjZ are preferable for 
the entanglement amplification is understood by considering the role of the measurement 
operator M + in the present protocol. When the measurement strength fulfills ([6]), M + 
reduces components with respect to |1) of an input state p (e.g., (ll|p|ll)). A large 
negative value of A SjZ implies that tr [(|1)(1| CgJ-^Ps] is a predominant component of the 
separable part of the input state. Therefore, the measurement operator M + considerably 
decreases the weight of the separable part in the input density matrix. We remark that 
positive values of A 8jZ are preferable for the entanglement amplification when a > (3 and 
the measurement strength satisfies e = a 2 — /3 2 (> 0). 

5. Implementing the protocol with general measurements 

We mention that no specific type of weak measurements is essential for this approach. 
We can use a generalized two-outcome measurement, 



with < p < 1 and < q < 1. We note that in the Procrustean method one of the two 
parameters is fixed as 1. The two real parameters p and q may be determined by the 
requirements \ip) = |^mse) and A + < A, given as, respectively, 



For pure initial states, the latter condition is not necessary. Adjusting the two real 
parameters depending on the initial partially-entangled state and the separable part, 
we may achieve entanglement amplification with some probability. Alternatively, we 
may use a continuous weak measurement. Continuously monitoring system A with very 
weak measurement strength, the entanglement of formation of the total system may 
fluctuate and take its maximal value 1 at some point in time. We note that both of 
them have the same efficiency as the Procrustean method. 

6. Summary 



Mi = Vp|o)(o| + ^|i)(i 



(27) 
(28) 




pa 2 = q (3 2 y 



(29) 
(30) 



(j>-q)(y/l-S(ip) + A,z)<0. 



We have proposed a local probabilistic protocol to make a maximally-entangled 
state from a partially-entangled pure state by applying elementary two-outcome weak 
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measurements. The present approach requires only the implementation of local 
measurements and the information of the initial entanglement. The probability to 
find the maximally-entangled state at each step decreases to zero with increasing 
measurement number. In addition, we found that the present approach corresponds 
to a decomposition of a two-outcome POVM which generates maximally-entangled 
or separable states as the measurement results. The success probability of the map 
corresponding to this POVM is high if the initial entanglement is large. The case for a 
mixed initial state was also examined. We found that the method for a pure state still 
leads to entanglement amplification, but not purifying a maximally-entangled state, if 
the separable part satisfies a certain condition. Several future issues will contribute to 
the development of weak measurement-based quantum control and the illustration of 
interesting and useful properties of POVM's. 

An interesting application example of this method is a quantum network system in 
which system A is a mediator (e.g., photon) and system B is a quantum memory. In 
this case, it is not easy to perform controlled operations between the two systems after 
their initial preparation. A local-operation-based approach will be effective for such a 
situation. 
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